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Abstract. Atoms can scatter light and they can also amplify it by stimulated 
emission. From this simple starting point, we examine the possibility of realizing 
a random laser in a cloud of laser-cooled atoms. The answer is not obvious as both 
processes (elastic scattering and stimulated emission) seem to exclude one another: 
pumping atoms to make them behave as an amplifier reduces drastically their 
scattering cross-section. However, we show that even the simplest atom model 
allows the efficient combination of gain and scattering. Moreover, supplementary 
degrees of freedom that atoms offer allow the use of several gain mechanisms, 
depending on the pumping scheme. We thus first study these different gain 
mechanisms and show experimentally that they can induce (standard) lasing. 
We then present how the constraint of combining scattering and gain can be 
quantified, which leads to an evaluation of the random laser threshold. The 
results are promising and we draw some prospects for a practical realization of a 
random laser with cold atoms. 
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1. Introduction: atoms as scatterers and amplifiers 

Among atom-light interactions, elastic - or Rayleigh - scattering is one of the simplest 
processes [l]. It has been known for a long time that in optically thick atomic 
vapors, multiple scattering leads to diffusion of light, or "radiation trapping" [21 [3]. 
However, only the progress in laser cooling and trapping [4] allowed the more recent 
demonstration and study of this effect with only true elastic scattering, by suppressing 
the Doppler- induced frequency redistribution fHIHllT]. Many experiments have also 
been performed to study the coherence properties of multiply-scattered light in such a 
medium, especially using coherent backscattering as a probe [8l[9] (see |10 [ [TT | [T2 t [T3] 
for review articles devoted to this topic). 

Cold atoms have also driven a renewed interest in non linear spectroscopy, as the 
suppression of Doppler broadening allows the resolution of narrow spectral resonances 
in pump-probe spectroscopic schemes [14]. One example is Raman gain between 
Zeeman sublevels [HI [16]. The idea of building a laser upon this gain followed soon 
after its first observation [IT]. Nevertheless, cold atoms in optical cavity have been 
then mainly used for quantum optics purposes [18] and laser demonstrations with 
different gain mechanisms are much more recent \19l I20j . 

If one could combine both gain and radiation trapping at the same time in a 
cold atom cloud, this could give rise to a diffusive random laser, as predicted by 
Letokhov [21]. Since his original "photonic bomb" prediction, great efforts have 
been made to experimentally demonstrate this effect in different kinds of systems 
[22 , 23, 24 , 25'j[26l[27l[28], as well as to understand the basic properties of random lasing 
[29. 30,, Jl, 32]. The broad interest of this topic is driven by potential applications 
(see [33] and references therein) and by its connection to the fascinating subject of 
Anderson localization [34l [35]. State-of-the-art random lasers [33l [Sfil \37\ [38] are 
usually based on condensed matter systems, and feedback is provided by a disordered 
scattering medium, while gain is provided by an active material lying in the host 
medium or inside the scatterers. In general, scattering and gain are related to different 
physical entities. 

The peculiarity of a random laser based on cold atoms would be that the same 
microscopic elements (the atoms) would provide both ingredients (scattering and gain) 
of random lasing. On the one hand, it leads to an easier characterization and modelling 
of the microscopic properties of the system, which can be extremely valuable for a 
better understanding of the physics of random lasers. Moreover, cold atoms are "clean" 
and perfectly characterized samples. In addition, relaxation of optical coherences is 
limited by radiative processes, which makes the transition to cooperative emission, 
such as superfiuorescence [39[ 140). more accessible than in condensed matter systems. 
On the other hand, it is clear that pumping atoms to induce gain reduces drastically 
their scattering cross-section, due to the saturation effect [l|: atoms spend less time 
in their ground state, in which they can scatter Hght. It is thus not obvious at all that 
reasonable conditions for random lasing can be obtained in cold atoms. Moreover, 
these conditions are expected to be different for each gain mechanism. 

The purpose of this article is to present the status of our experimental and 
theoretical investigations on this issue. In the next section we present our work 
on lasing with cold atoms, with a standard cavity [20]. Experimentally, this 
demonstration of a cold-atom-based laser with different gain mechanisms is a first 
important step towards the building of a random laser. More fundamentally, this 
experiment can be compared to theory to validate or improve the modelling of such 
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gain media, including saturation effect, laser dynamics or statistics. We present a part 
of this modelling in section [31 Ultimately, when studying the random laser properties, 
it will be of first importance to be able to discriminate which behaviour originates 
from the particular used gain medium and mechanism and which one is specific to the 
feedback mechanism (standard cavity or scattering) . The cavity-laser thus serves as a 
reference, to be compared with theory and with the (forthcoming) random laser. 

The second important step of this project, presented in section IH is to quantify 
the constraint of combining gain and scattering at the same time. This leads to the 
evaluation of the random laser threshold |4H I42j . The goal is first to establish the 
feasibility of a random laser with cold atoms, which was not obvious, and second 
to compare the different gain mechanisms in order to experimentally choose the 
more appropriate one with the best pumping parameters. This is necessary before 
conducting further experimental efforts to achieve the threshold. As the results are 
promising, we draw some prospects for our future work in the last section. 

2. Gain and lasing with cold atoms: experimental investigation 

In this section, we briefiy describe our experimental setup and then the different 
gain mechanisms we have studied. Mollow gain is the simplest as it involves only 
one pumping field and two- level atoms O |44] . By using the more complex atomic 
structure of rubidium atoms, we can create two-photon transitions between two non- 
degenerate ground states. This can produce Raman gain \TE \ fTB t flT]. Finally, the 
atomic non-linearity can give rise to parametric gain, for example with four-wave 
mixing. We have demonstrated laser action with each of these mechanisms (figure [T| 
[20j . We mention also briefiy some other gain mechanisms that could be used. 
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Figure 1. Transverse modes of cold-atom lasers, (a) Gaussian TEMqo mode, 
obtained by inserting a small diaphragm in the cavity. Typical modes of: (b) the 
Mollow laser, (c) the Raman laser, and (d) the four-wave mixing laser respectively. 

2.1. Experimental setup 

Our experiment uses a cloud of cold ®^Rb atoms confined in a vapour-loaded Magneto- 
Optical Trap (MOT) [4] produced by six large independent trapping beams, allowing 
the trapping of up to 10^° atoms at a density of 10^" atoms/cm'^, corresponding 
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to an on-resonance optical thickness 60 ~ 10. A linear cavity, formed by two mirrors 
separated by a distance L — 0.8 m, is placed outside the vacuum chamber. Reflections 
on the vacuum cell yields a low flnesse = 16. To add gain to our system, we 
use either one or two counter-propagating pump beams, denoted F (forward) and B 
(backward), produced from the same laser with a waist Wpump — 2.6 mm, with linear 
parallel polarizations and a total available power P = 80 mW, corresponding to a 
maximum pump intensity / = 2P/{TTWp^^^p) « 750 mW/cm^. The pump is tuned 
near the F — 3 ^ F' — A cycling transition of the D2 line of ^^Rb (frequency ujq, 
wavelength A = 780 nm, natural linewidth r/27r = 6.1 MHz), with an adjustable 
detuning A — o^f.b — ^^o and has an incident angle of about 20° with the cavity axis. 
An additional beam P is used as a local oscillator to monitor the laser spectrum or 
as a weak probe to measure transmission or reflection spectra (flgures IDJb) andlUJb)) 
with a propagation axis making an angle with the cavity axis smaller than 10°. Its 
frequency up can be swept around the pump frequency with a detuning S — wp — wf,b- 
Both lasers, pump and probe, are obtained by injection-locking of a common master 
laser, which allows to resolve narrow spectral features. In our experiments, we load a 
MOT for 29 ms, and then switch off the trapping beams and magnetic fleld gradient 
during 1 ms, when lasing or pump-probe spectroscopy are performed. In order to 
avoid optical pumping into the dark hyperflne F = 2 ground state, a repumping laser 
is kept on all time. 

2.2. Gain mechanisms 

Contrary to most laser gain media, cold atoms do not present non-radiative, fast- 
decaying transitions, preventing a standard 4- level scheme to produce a population 
inversion, as e.g. in many gas lasers. Nevertheless there are many different mechanisms 
which allow an inversion between two different atomic states. These can either be 
different states in external degrees of freedom (momentum or vibrational levels in an 
external potential) or internal degrees of freedom (dressed states [T] or different ground 
states). We have used two mechanisms, which we call respectively Mollow and Raman 
gain (described in the two next subsections), to induce lasing with a cavity and to 
study the threshold of a random laser (section [4. 3p . The atomic nonlinearity can also 
be used to obtain parametric gain and lasing (subsection 12. 2. 3p . as well as other, more 
complicated, schemes using quantum interferences or the external degrees of freedom, 
as briefly mentioned in subsection 12.2.41 

2.2.1. Mollow gain The most simple gain mechanism we can imagine in cold atoms 
was described by Mollow [43] and observed soon afterwards [44] with atomic beams. It 
involves a two-level atom driven by one strong pumping field. The driving field induces 
a population inversion in the dressed-state basis [l] (see figure [7l[ a)) and therefore a 
weak probe beam can then by amplified. The whole process can also be described in 
the bare-states basis by a three-photon transition from the ground state to the excited 
state via two absorptions of pump photons as sketched in figure [2]Ja) . 

The main amplification feature appears for a pump-probe detuning 6 = 
sign(A) vA^Pf^, where Q, is the Rabi frequency of the pump-atom coupling, related 
to the pump intensity / by il^ = C^r^//(2/sat) (-^sat = 1-6 mW/cm^ is the saturation 
intensity and C is the averaged squared Clebsch-Gordan coefficient of the F = 3 ^ 
F' = A transition) and has a typical width on the order of F. Note that another, 
dispersion- like feature appears around (5 = 0, which is associated with two-photon 
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Figure 2. (a) Principle of the MoUow gain depicted as a three- photon transition 
from the ground state to the excited state. It can also be viewed as a population 
inversion in the dressed-state basis (figure Efa)). (b) Transmission spectrum, 
computed for typical experimental parameters bo = lOi f2 = SF and A = F. 



spontaneous emission processes [45] . This contribution also induces gain but with 
a much smaller amplitude. Note that this can generate lasing without population 
inversion [46l[47l[48] . 

In our experiment, we have measured single-pass gain as high as 50 %, which is 
more than enough to induce lasing even with a low-finesse cavity [20j. This Mollow 
laser has an output intensity reaching 35 /iW achieved for |A| ~ 2r. Its threshold 
in pump intensity is in agreement with the corresponding measured single-pass gain 
and the losses of the cavity. The laser polarization is linear, parallel to the pump 
polarization, because this is the configuration for which gain is maximum, as the 
driven atomic dipole is then parallel to the probe field. We have also measured the 
output power as a function of the pump intensity, as reported in figure [3l We observe 
a threshold and the maximum intensity is reached for a pump power about 5 times 
larger than at the threshold. This behaviour is well described by the theoretical model 
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Figure 3. Laser power versus pump power. Points are experimental data 
obtained after an average of 1000 experiments, the continuous line is the 
theoretical fit (see section [Sj with a single fit parameter rj ~ 8 (corresponding 
to iV ~ 10* atoms in the cavity mode) and all other parameters fixed at their 
experimental values, 6o = H, A = F, = 16, L = 0.8 m, and cavity waist 
at the MOT location w = 500 /im. The horizontal axis is normalized by the 
threshold value of the squared Rabi frequency fi^j^j, ; the vertical axis is normalized 
by the maximum value of the electromagnetic field intensity reached at about 

(n/Qthr)' ~ 5. 
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presented in section [31 

The Mollow laser works for pump detuning |A| < 4r. When the pump frequency 
in detuned farther away from the atomic resonance, Raman gain becomes dominant, 
and the system switches to another regime of laser, based on Raman gain. 

2.2.2. Raman gain Raman transitions refer in general to two-photons transitions 
between two non-degenerate ground states, the intermediate energy level being in the 
vicinity of the atomic excited states. To obtain gain, a pumping field can induce the 
first upward transition and a probe beam can then be amplified at the frequency of the 
downward transition. In our case, Raman gain relies on the pump-induced population 
inversion among the different light-shifted mp Zeeman sublevels of the F — 3 hyperfine 
level [TSl [TB] , as depicted in figure SJa) . The optical pumping induced by the tt- 
polarized pump laser leads to a symmetric distribution of population with respect to 
the mp = sublevel of the ground state, with this sublevel being the most populated 
and also the most shifted, due to a larger Clebsch-Gordan coefficient f49|. To record 
a transmission spectrum, atoms are probed with a 7r-polarized (with perpendicular 
direction) probe beam, thus inducing Amp = ±1 Raman transitions. Depending 
on the sign of the pump-probe detuning S, the population imbalance induces gain 
or absorption. Each pair of neighboring sublevels contributes with a relative weight 
depending on the population inversion. In practice however, the contributions of 
different pairs are not resolved and only two structures (with opposite signs) are 
visible, one corresponding to amplification for S — — (5r and one to absorption for 
(5 = i5r. Note that this situation corresponds to a red detuning for the pump (A < 0) 
and that the signs are inverted for blue-detuning (A > 0). As Sn comes from a 
differential Hght-shift (because of different Clebsch-Gordan coefficients), it is usually 
on the order of F/IO, whereas A is a few T. The width 7 of the resonances is related 
to the elastic scattering rate, also much smaller than F [161 . Far from the main atomic 
absorption resonance, the Raman resonance is thus a narrow spectral feature, as in 
figure 

The laser obtained with Raman gain has an output polarization orthogonal to 
the pump one (contrary to the Mollow laser) and less power (2 /iW). Moreover, the 




Figure 4. (a) Principle of the Raman mechanism, depicted here for a F = 
1 — > -F' = 2 transition, (b) Experimental transmission spectra, plotted as a 
function of the pump-probe detuning <5. Without pumping, spectrum (1) shows 
only the atomic absorption. A pump beam of detuning A = —3.8 F and intensity 
13 mW/cm^, corresponding to a Rabi frequency f! = 2.5 F, is added to obtain 
spectrum (2), which then exhibits a Raman resonance in the vicinity of (5 = 0. 
Moreover, the atomic absorption is shifted due to the pump-induced light shift 
and the absorption is reduced due to saturation. 
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sharpness of the gain curve makes the Raman gain very sensible to any Doppler shift. 
The radiation pressure from the pump beam makes thus the laser emission to stop 
after only ^ 20 ^s. On the other hand, the narrow spectrum of the laser can be easily 
characterized by a beat-note experiment |20j . 

2.2.3. Parametric gain The Mollow and Raman lasers only require one pump 
beam. By using a second pump beam, we can induce four wave mixing (FWM): 
the two pumps of frequencies ujf and cob and one probe - or initial fluctuation - 
of frequency ujp generate a fourth field at frequency ojq, called the conjugate field 
|50[ [5T| I52j . The frequencies and wave- vectors of all the fields are related by energy 
and momentum conservations. If we want to obtain gain for the probe, we have to 
choose a configuration where the conjugate frequency equals the probe one: ujc — top. 
Then, the pump frequencies have to fulfill the condition cup + lub = 2u!p. From an 
experimental point of view, the most simple configuration consists of all frequencies 
to be the same ("degenerate FWM"). This is the experimental situation that we have 
studied so far, and we did obtain lasing in cold atoms with this mechanism [2^. Note 
that this mechanism has been observed a long time ago with hot atoms [53t [54j [55] . 

Due to the phase matching condition, however, the gain is not in the forward 
transmission of the probe beam, but in backward refiection, provided that the two 
pumps are counterpropagating. The conjugate beam is actually the phase-conjugate 
of the probe beam. This property has a number of consequences for the laser [20] . 
First, it leads to a different threshold condition [54]: a refiectivity of only 1 % is enough 
to generate lasing, despite the much larger losses of the cavity (32 % for a round trip). 
This is due to constructive interferences between transmitted and refiected waves, as 
observed in double-pass experiments [56j. Second, it leads to more complex transverse 
modes (figure [T]Jd)), because the phase conjugation mechanism allows any transverse 
pattern to be stable through the resonator [57j . 

Finally, these properties lead to a much larger power than with Mollow and 
Raman gain, as up to 300 have been obtained. Moreover, this laser might find 
appHcation for other topics Hke pattern formation [58l [59] or the production of twin 
beams for quantum optics [60l [6^ and even quantum imaging [62j. Note also that 
the use of two pump beams allows a longer laser emission (figure [6]Jb)) because the 
effect of radiation pressure on the cloud is suppressed (other effects can arise due to 
the dipole force, see [63j). 




Figure 5. (a) Principle of four-wave mixing, (b) Typical experimental reflection 
spectrum. 
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2.2.4- Other gain mechanisms Our previous study on the gain mechanisms that can 
be used in cold atoms is not exhaustive. For example, the two hyperfine ground states 
of rubidium atoms can also be used to produce Raman gain [64] . 

Other, more complicated, schemes, involve quantum interference to induce gain 
without population inversion (whatever the basis) [48]. This can be realized with a 
A scheme [65] or a V scheme [Ml [67]. In this last configuration, a large detuning 
between the pump and the gain frequency can be reached by using the two D lines of 
rubidium, whose separation is 15 nm. If used to produce a random laser, this can be 
highly valuable to make the detection easier (see the discussion of section [U . 

Another possibility is to use the atomic external degrees of freedom, i.e., their 
kinetic energy. Transitions between different velocity classes produce recoil-induced 
resonances [68], and high gain can be achieved [69]. These resonances can ultimately 
lead to a "Collective atomic recoil laser" [TQl [71] , which has been demonstrated with 
cold atoms |19) . 

Finally, one could also consider higher-order photonic processes, such as two- 
photon dressed-state lasers [72] . 



2.3. Laser dynamics 

We have been interested so far only in steady-state properties. Our cold-atom-based 
laser could also reveal interesting dynamical properties. 

For example, we have observed that the Mollow laser switches on very quickly, 
whereas the FWM laser switches on very slowly (figure [6] - notice the change in the 
time scales). The laser rise-times are related to the spectral properties of the gain 
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Figure 6. Temporal evolution of the Mollow (a) and FWM lasers (b). The 
pump beams are switched on at t = 0. Note the change in the time scales. The 
Mollow laser switches on very fast (the initial decrease is due to optical pumping 
to the dark hyperfine state), whereas the FWM laser needs more time. Inset of 
(b): Power spectrum, with a logarithmic scale, of the FWM laser measured by 
a beating with a local oscillator of frequency uq. All data are the result of an 
average of 1000 cycles. 
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mechanism. We are indeed in the bad cavity hmit, where the decay rate of the cavity 
is not smaller than the spectral width of the amplification. The measured orders of 
magnitude (~ 0.2 for the Mollow laser, more than 100 /iS for the FWM laser) are 
consistant with the inverse of the laser spectral width (a few MHz for the Mollow 
laser, a few kHz for the FWM laser). Note that the spectral width of the FWM 
laser, measured by a beating with a local oscillator (inset of figure E^b) - note the 
logarithmic scale) , is substantially smaller than the width of the reflectivity spectrum 
(figure mb)). 

Clearly, further experimental and theoretical studies would be needed to fully 
understand the dynamical behaviour of these different lasers. However, a first insight 
is given in the next section. 



3. Model of the Mollow laser 



We present in this section a theoretical analysis of our cold-atom laser based on Mollow 
gain. 



3.1. Formalism 



We consider an ensemble of N identical two-level atoms pumped by a strong field of 
frequency ujl, detuning A and intensity characterized by the Rabi frequency fl. They 
are placed inside a cavity of eigenfrequencies uJm, and Um, aJn are the annihilation and 
creation operators, corresponding to the mode m. The cavity is partly open and the 
energy of the electromagnetic field can escape at a rate 2k. We recall that we are in 
in the bad cavity limit [73] . 

To describe the dynamics of the coupled system of atoms and modes of the 
electromagnetic field, we use the master equation approach [Tl[73l[74]. The density 
matrix operator obeys 

p = -i[H,p]+ LaP + LfP , (1) 
where we set ft = 1 and the Hamiltonian is 



H = 



1 ^ 



-A 



n 



M 



(ffimcrja™ + h.. 



M 

E 

m— 1 



+ 6ma\rfim ■ (2) 



Here the rotating- wave approximation was used and the operators a^j, aj and aj 
are standard Pauli matrices, describing the populations {(Tsj) and the coherences (cj) 
of the j-th atom. The detuning of the mode frequency from the pump frequency 
is 6m and gjm are the coupling constants between atoms and cavity modes, which 
we assume to be independent of the mode index m and to have the same absolute 
value g for all atoms and random phases 4>j: gjm — 5'exp(i(/)j). In the Hamiltonian 
of equation only a part of modes of the electromagnetic field in the cavity {M 
modes) are treated explicitly, whereas spontaneous emission into other modes is taken 
into account through the first non-Hamiltonian term of the master equation |(T]) : 



F ^ 



{2a,pa] 



paUj 



(3) 



Towards a random laser with cold atoms 



10 



The second non-Hamiltonian term in equation ([T]) describes cavity damping of the 
electromagnetic modes: 

M 

Lpp = «^ X! i'^^'rnPaln ^ al^amP - Pa-lnO-m) ■ (4) 
m— 1 

The physical processes described by the set of equations become especially 

clear in the dressed-state basis [I]. The dressed states are eigenstates of the system 
"atoms + pump field". In the dressed-state basis, the states |1)^ and |0)^ of the 
two-level atom j become |-t-)^ = (cos0) + (sin6') |0)^- and = —{sind) + 

(cos 6*) |0)^., respectively, where Q = f7'sin26' and A = - ft' cos 29 with = ^/WT^ 
being the generalized Rabi frequency (figure [Tlja)). Equations ^ and ^ can then be 
transformed (equation |[4]) remains unchanged) and equation |(T]) leads to the following 
semi-classical equations for quantum- mechanical expectation values of operators aj , 
CT3j and [731174]: 

CT3j = - 7i {(^3j - 0-3j) + *(crj, a„) , (5) 
o-j = - (72 + i^') CTj + ^'(crj, (Tsj, a„i) , (6) 
flm = - {K + i5m)am + Q{(J3j,o-j) , (7) 

where asj = -2cos26'/(l-hcos2 261), 71 = (r/2) (l + cos ^ 26>), 72 = (r /4) (2 + sin^ 261) 
and the expressions of functions $, ^, are given in [Appendix A[ The equations 
for (7* and are complex conjugates of equations ^ and respectively. These 
equations describe respectively the evolution of the populations, the coherences and 
the cavity field. The first term of each right-hand-side contains the natural evolution 
associated with rates 71, 72 and k as well as the driving by the pump field (Rabi 
frequency fl). The equations are coupled by their second r.h.s. terms (see equations 
lA.lHA.Sl for the complete expressions), which contain the atom-field coupling g and the 
number N of atoms coupled to the cavity mode. For a Gaussian cavity mode with waist 
w one obtains [75] g ~ y^ScT/L/kow, N ~ {kQw)^bo/12 where fco — ujo/c — 2n/X, 
bo is the optical thickness of the cloud at resonance without pump and L the length 
of the cavity. To take into account the fact that volumes and shapes of the atomic 
cloud and cavity modes are different and not Gaussian, we will multiply by a free 
parameter rj that will be adjusted to fit experimental results. 

Equations ([SHll) can be solved numerically following the procedure described in 
[Appendix A 

3.2. Dynamics of the laser 

The dynamics of the laser depends strongly on how the different time scales associated 
with the evolutions of populations, coherences and field can be compared. Following 
the classification introduced in |76[ I77j. class A lasers are those in which population 
and polarization decay much faster than the field so that the dynamics is governed by 
the single field equation (for example He-Ne, Ar+, Kr+, dye lasers), class B lasers are 
those for which the population decays slowly and two equations are necessary (e.g. 
ruby, Nd, CO2 lasers) and finally, class C lasers are those for which the three time 
scales are of the same order of magnitude. The last ones are known to exhibit chaotic 
behaviour [761177] . 

For our cold atom sample, the relaxation rate of atomic coherences is only given 
by optical processes (no coUision, e.g.). As a consequence, and with Mollow gain. 
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Figure 7. (a) Coupling of the bare atomic states \0,n), |l,n) (left) to the dressed 
states \ — ,n), \ + ,n) (right), where n stands for the number of pump photons. The 
population inversion and the subsequent gain is sketched on the dressed-state 
representation, (b) Laser power versus pump power for a laser with a single 
cavity mode interacting with atoms (M = 1). Values of parameters are chosen 
as in the experiment, bg = 11, A = F, = 16, L = 0.8 m, ui = 500 /im, and 
we have taken Sm = l.SF and rj ~ 15, corresponding to N = 2.3 x 10®. Regions 
I, II and III correspond to three different types of behaviour illustrated in figure 
[3 Regions III (chaotic behaviour) appear only when rj is large enough, (e.g., no 
chaotic behaviour is observed for r; = 8 as in figure [3j. 



the time scales for the evolution of the populations and coherences are both on the 
order of (note that this is not the case with Raman gain or FWM). The cavity 
damping rate k is related to the finesse J- of the cavity by k = Ticj^LT . With our 
experimental parameters, it gives K/27r ~ 5.9 MHz, i.e., k ~ F and we are thus 
in a regime corresponding to class C lasers. This is consistent with the behaviour 
observed in our numerical simulation, as described in figures [7l[b) and [HI We show 
the dependence of the laser power on the pump power in figure (T^b) for the case of 
when only a single mode (M = 1) with a well-defined frequency (detuning 6 — l.SF) 
is present. The non-monotonic character of this dependence is accompanied by a 
complex dynamic behaviour that we illustrate in figure [8l Three distinct dynamic 
regimes are possible, depending on the pump strength: the laser emits at a single 
frequency (I) or at several frequencies at a time (II). In addition, chaotic behaviour 
appears for certain (HI)- These results show that Mollow gain, associated to a low- 
finesse cavity, is perfectly suited for the study of chaotic regimes in lasers, as already 
predicted in [73]. 

We did not observe this kind of dependence of laser power as a function of the 
pump intensity in our experiment. This behaviour appears indeed only for some values 
of atom number or detuning (S^ and a specific experimental implementation would be 
needed. To model the experimental results, we solved the coupled laser equations 
with a well-defined single mode of the field (M = 1, detuning &\ — S) and average 
the resulting laser intensity |ap over 200 equidistant (5 between — 2F and SF. This is 
supposed to account for the fact that the exact frequency of the cavity mode is not 
precisely known in the experiment and fiuctuates shot to shot, as well as for averaging 
over many consecutive measurements (1000 in the experiment). Using 77 as a free 
parameter, we are then able to obtain a good fit to experimental data (solid line in 
figure [3]). The obtained value of 77 = 8 is in qualitative agreement with the observation 
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rt > 100. The parameters used in the calculation are the same as in figure [T] 



of high-order transverse laser modes: the corresponding cavity modes have a larger 
overlap with the atom cloud compared to a gaussian TEMqo mode. 
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4. Threshold of a random laser with cold atoms 

We now turn to the more specific question of realizing a random laser based upon 
one of the gain mechanisms studied in the previous section. As explained in the 
introduction, gain is not the only necessary ingredient for building a random laser, as 
scattering should also be present to provide feedback. In the previous laser experiment, 
scattering was not considered, and the optimum conditions for the laser operation are 
those which optimize the single-pass transmission though the cloud, which are very 
probably not the one that optimize scattering. To establish the best conditions for 
random lasing, we thus have to quantify more precisely the constraint of combining 
gain and scattering. This condition is already contained in the seminal Letokhov's 
paper [21], and leads indeed to the evaluation of the random laser threshold. 

4.1- Random laser threshold 

Although different kinds of random lasers have been achieved experimentally, several 
points are still under debate. Random lasers may be divided into two categories, 
depending on the feedback mechanism [37j. In the regime of incoherent (or intensity) 
transport, feedback is provided by an increase of the photon path lengths (or lifetime) 
in the system. This regime has been observed experimentally in many systems, such 
as powdered laser crystals [22], dye solutions containing microparticles |23j, soHd- 
states materials [38], or porous glass infiltrated by a dye solution in liquid crystal [26j . 
In a regime in which interference effects survive the random scattering process, a 
coherent (or field) feedback is expected, leading to a behaviour closer to that taking 
place in conventional lasers. This behaviour is expected to become substantial close 
to the Anderson localization regime [34]. Experiments in this direction have been 
carried out, e.g., on zinc oxide powders [MllS^, and signatures of coherent feedback 
on the random lasing mechanism have been discussed recently [3T| [32} [35] . For a given 
practical realization of a random laser, discriminating between a coherent or incoherent 
feedback mechanism is a complicated task. For example, narrow emission peaks above 
threshold can be observed, even far from the localization threshold [37l [78]. The 
existence of such peaks in experiments can be attributed to exponential gain along 
very long diffusion paths [78], or to interferences and coherent feedback in the weak 
scattering limit [31]. On the route towards a random laser with cold atoms, we first 
consider a diffusive random laser (with intensity feedback) . This approach is motivated 
by the fact that experiments that should be carried out in the near future are expected 
to work in this regime. Another motivation is the identification of specific signatures 
of the incoherent feedback regime that could guide the experimental observations. A 
deviation from these signatures could be a measure of the onset of a coherent feedback 
mechanism. 

4-2. Letokhov's threshold for cold atoms 

From Letokhov's diffusive description of fight transport in a homogeneous, disordered 
and active medium of size L, we know that the random laser threshold is governed 
by two characteristic lengths: the elastic scattering mean free patlQ -^sc [79] and the 
linear gain length £g {£g < corresponds to absorption or inelastic scattering). In 

* We consider only isotropic scattering so that the transport length equals the scattering mean free 
path. 
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the diffusive regime, defined as L ^ 4c, the lasing threshold is reached when the 
unfolded path length, on the order of L^/4c, becomes larger than the gain length. 
More precisely, the threshold is given by [211 IM] 



LoS > Pn^feZiJi , (8) 

where /? is a numerical factor that depends on the geometry of the sample (/3 = 1 for a 
slab, /? = 2 for a sphere, which is the case we consider in the following), and Left = rjL 
is the effective length of the sample, taking into account the extrapolation length |79j . 
For L > 4c and a sphere geometry, 77 = 1 + 2^/ [L/4c + 2^] with ^ - 0.71 [8OI [81]. 
Note that deeply in the diffusive regime (L ^ 4c), ?7 1- Another important length 
scale is the extinction length 4x, as measured by the forward transmission of a beam 
through the sample, T — e'^^^"". The extinction length is related to the other lengths 
by = l^^ ~ ■ Note that this reasoning may not be appropriate for backward 
gain such as produced by FWM. 

For an atomic vapour, these characteristic lengths can both be computed as a 
function of the atomic polarizability a{uj) at frequency w. The extinction cross-section 
is indeed given by tTex(w) = k x Im[a(a;)] and the elastic scattering cross-section by 
CTsc(w) — k^/dw X |a(aj)p [82] {k = lojc is the wave vector). Note that the first 
relation is general to any dielectric medium whereas the second one is specific to dipole 
scatterers. The characteristic lengths are then ^^^sc — po'cx,sc, where p is the atomic 
density. The gain cross-section can be defined the same way by ^ — pag. The vapour 
is supposed at constant density and homogenously pumped, so that both p and a are 
position-independent. Even though this is not the precise geometry of a cold-atom 
experiment, it allows us to perform analytical estimations. As we consider resonant 
scatterers, we deal only with quasi-resonant light and we shall use k = ko ^ uq/c 
with luq the atomic eigenfrequency. In the following, we shall also use a dimensionless 
atomic polarizability a, defined as a = a x 67r/fcQ, and omit the dependence on uj. 
We can now rewrite CTsc = cro|<5p and Cg = (Tq (jap — Im(a)), where (Tq — 6n/kQ is 
the resonant scattering cross-section (for aJ = O^J=l transition), such that the 
threshold condition, as expressed by equation ([8]), reduces to [41] 

27r 

paoLcs = ?70o > — = , (9) 

V3|a|2 (|5|2-Ini(5)) ^ ' 

where bo is the on-resonance optical thickness of the cloud. This condition is valid as 
soon as the medium exhibits gain, i.e., jap — Im(Q;) > 0. Interestingly, the condition 
Im(Q!) < 0, corresponding to single-pass amplification (T > 1), is not a necessary 
condition. 

The threshold condition is thus given by a critical on-resonance optical thickness, 
which is an intrinsic parameter of the cloud, expressed as a function of the complex 
atomic polarizability only, which depends on the pumping parameters. Although the 
initial condition of equation ^ involves two characteristic lengths, we emphasize here 
that this is really one single independent parameter, as real and imaginary parts of 
the atomic polarizability are related via Kramers-Kronig relations ^83j. This point 
is due to the originality of the system that we are considering, in which the same 
atoms are used to amplify and scatter light. This property can be fruitfully used to 
experimentally determine the threshold, as only one single measurement can provide 
enough information. A weak probe transmission spectrum, which we can rewrite with 
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our notations, 

T(cj) e-'"'^™I"('^)l , (10) 
allows indeed the full characterization of ce{uj). 

4-3. Application to Mollow and Raman gain 

We now apply the previous results to Mollow gain and to Raman gain. For the first 
one we use the ab initio knowledge of the polarizability of a strongly pumped two-level 
atoms. For the second one, we use experimental transmission spectra T{u!) to extract 
the atomic polarizability via equation (fTO|) and Kramers-Kronig relations. 

With Mollow gain, the polarizability is exactly and analytically known [43], 
assuming a weak probe field, which is a good hypothesis for calculating the threshold 
of the random laser: 



1 + 4A2 



2. /.OA .^ (11) 



2 1 + 4A2 + 2172 

{S + t){S - A + i/2) - n^S/{2A - i) 
{6 + i){S - A + i/2){S + A + i/2) ^n^{S + i/2) ' 



where A, S and Q, are in unit of F. 

For each pair of pumping parameters {A,il}, the use of the polarizability (fTTj) 
into the threshold condition ^ allows the calculation of the critical on-resonance 
optical thickness bo as a function of the pump-probe detuning S. Then, the minimum 
of 6o and the corresponding 6 determine the optical thickness 6ocr that the cloud must 
overcome to allow lasing, and the frequency (Srl of the random laser at threshold. The 
result is presented in figure [DJa) for a spherical geometry {(3 = 2). The result for 6ocr 
is independent of the sign of A. The minimum optical thickness that allows lasing is 
found to be 6ocr ~ 200 and is obtained for a large range of parameters, approximately 
along the line fl « 3 A. The optimum laser-pump detuning is near the gain line of 
the transmission spectrum, i.e., (5rl ~ sign(A) vA^PfP, with however a small shift 
compared to the maximum gain condition due to the additional constraint of combined 
gain and scattering. 

With Raman gain, we do not have a simple analytical expression of the 
polarizability. Nevertheless, we can measure its imaginary part by a transmission 
spectrum (equation (fTO|) ). A fit of the result, followed by its transformation by 
Kramers-Kronig relations [83], allows us to recover the full complex polarizability 
and to use equation ^ [42] . Note that this method is general and could be appHed to 
any gain mechanism. However, Raman gain is perfectly suited for this analysis as we 
can fit the corresponding transmission profile by two inverted Lorentzian lines (see the 
Raman structure on figure|4]Jb), near 6 = 0), the transformation of which via Kramers- 
Kronig relations is analytic and well-known (it gives a standard dispersion profile) . As 
the two Lorentzians are not well separated (for most parameters, (5r < 7), it leads for 
the scattering cross-section (cx |a|2) to a bell-shaped curve centered near 5^0 [42\. 
Note that besides creating gain, the Raman transition adds some scattering. 

As for Mollow gain, we have varied the pump parameters {A, $7} and for each pair, 
we have recorded a transmission spectrum and deduced the corresponding random 
laser threshold. We report the result on figure [9]Jb). The calculation contains 
corrections due to the peculiar polarization used in the experiment but deals only 
with the Raman resonance. It takes into account neither the main atomic absorption 
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(a) (b) 




Figure 9. Threshold of random lasing based on MoUow gain (left) and Raman 
gain (right) for each pair of pumping parameters A (detuning) and fl (Rabi 
frequency). We use the theoretical knowledge of the atomic polarizability Ijllll 
for MoUow gain whereas the prediction for Raman gain is based on experimental 
spectra (see text). In both cases, we obtain a minimum optical thickness of about 
200. 



line at ojq (which adds elastic scattering), nor the losses due to inelastic scattering. 
Doing so, the two corrections almost compensate each other, and leads to 6ocr ^ 200 
as the correct order of magnitude, obtained forO = A~3 — 4r j42] . 

These results are very promising. Indeed, by using well-established techniques to 
compress magneto-optical traps (see, e.g., [84^ I85j). achieving an optical thickness of 
200 should not be too difficult. 

4.4- Beyond the diffusive model 

The calculation of the lasing threshold based on Lethokov's approach leads, for the 
optimum pumping parameters, to ratio L/£sc = &ocr|fip ~ 0.5 for Mollow gain and 
L/isc ~ 2 for Raman gain. Because these values are not consistent a priori with 
the domain of validity of the diffusion approximation [S^, a more refined transport 
model is needed in order to check the relevance of the approach. A model for random 
lasers with incoherent feedback, going beyond the diffusion model, was introduced by 
Noginov et al. [87] , based on a phenomenological one-dimensional modelling of light 
transport. In this work, we use a more advanced three-dimensional model based on 
the Radiative Transfer Equation (RTE), that was introduced as a tool to predict the 
threshold of "classical" random lasers [88], and that we extended recently to describe 
a cold- atom random laser with Mollow gain [4T] . 

The RTE is a Boltzmann-type transport equation [89], that has a wider range 
of validity with respect to the ratio L/£sc than the diffusion equation [90]. The basic 
quantity is the specific intensity /j^(r, u, i), which describes the number of photons 
at frequency uj, at point r, propagating along direction u at time t. In order to 
compare the diffusive and RTE predictions, we will consider a slab geometry, with 
the z axis normal to the slab surfaces. In the case of plane-wave illumination at 
normal incidence, and for isotropic scatterers, the specific intensity only depends on 
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the space variable z and the angular variable fi — cos 9, with 9 the angle between the 
propagation direction u and the z-axis. In a system exhibiting gain and (isotropic) 
scattering, the RTE reads: 

- -^{z, ^l,t) + ^i -£(z, M, t) O m, t) + (2 4c)"^ j U^, m', t) d^l' 

(12) 

where c is the energy velocity in the medium. 

A feature of the RTE is that a modal expansion is available, whose asymptotic 
behaviour at large length and time scales leads to the modal expansion of the diffusion 
equation [80]. Therefore, under the conditions of uniform (in space) and constant (in 
time) gain, it is possible to build a modal theory of random lasers with incoherent 
feedback based on the RTE, that generaHzes Lethokov's approach beyond the diffusive 
regime [88]. We focus on the slab geometry (/? — 1) since the modal expansion of 
the RTE is well known in this case [80] (to our knowledge, no simple expansion is 
available for a sphere in the RTE approach) . The modal approach consists in looking 
for solutions of the form I{z,fi,t) — /k,s(m) exp(jKz) exp(si), where I{z,fi,t) is the 
specific intensity and the dependence on cu has been omitted for simpHcity. k can be 
chosen as a real parameter, and s can take complex values. For a given real k, s(k) and 
Ik,,s form a set of eigenvalues and eigenfunctions of the RTE. For isotropic scattering, 
eigenvalues and eigenfunctions can be obtained analytically [80|. In a passive medium, 
the eigenvalue corresponding to the mode with the longest lifetime read^: 

so{k)/c = £g ^ - - K/tan(K4c)] , for k4c < | • (13) 

In the presence of gain, the lasing threshold is reached when this modes starts to 
display an exponential amplification in time. Determining the threshold parameters 
amounts to calculating the gain length and scattering mean free path generating an 
eigenvalue so{k) > 0. Exactly at threshold, one has so{k) = 0. For a slab of width 
L, the dominant mode corresponds to k = n/LcS = Tr/{L + 2^4c)- In practice, the 
determination of k is meaningful as long as ^ = 0.71 can be taken as a constant 
(independent on L). Although not shown for brevity, we have verified with a full 
numerical solution of the RTE that this is the case as soon as L > 4c- This condition 
sets a limit of accuracy of the modal approach. 

It is well known that the diffusion approximation is asymptotically reached from 
the RTE in the limit of long time and large length scales. Indeed, in the limit k4c ^ 1, 
a first-order expansion of equation (fTSl yields 

s^^^\K)/c = £^'-K^esc/3, (14) 

which is the dispersion relation of the modes of the diffusion equation (DA) [80l lOT] . 
The associated threshold condition s^q^'^\k = tt/Lcs) = corresponds exactly to the 
result obtained with Letokhov's approach, with (3 = 1. A careful look at equations (flSl) 
and (fT4|) shows that the predictions of the RTE and diffusion approximation should 
not be substantially different. Firstly, we note that the gain contribution to so{k) 
is the same in both models (the first term is the same in equations (fT3|) and lfT4|) ). 
Secondly, the scattering contribution (the second term in the equations) is larger in 

tt Note that a misprint occurred in the published version of [41] and that the correct expression is 
given here. 
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Figure 10. Continuous line: Emitted random laser power normalized to the 
pump fluorescence power, as a function of the pump intensity. Dashed line: 
Normalized laser detuning <5r,l- The random medium is a spherical cloud of 
two level atoms with an on-resonance optical-thickness 6o = 600. 



the RTE model, but by a factor that remains smaller than 1.13 (when L ^ £sc)- 
This means that the correction introduced by the RTE remains relatively small, at 
least for the slab geometry. Numerically, it corresponds to an increase of ry^ocr of 
at most a few percents [41]. We can conclude that the model based on the diffusion 
approximation gives accurate results, even in the regime L ^ igc- Also note that the 
accuracy should be better with Raman gain, for which the optimum conditions are 
reached for L/£sc ''-^ 2 (instead of L/4c 0.5 for Mollow gain). Finally, we stress that 
our prediction shows that in a cold-atom system, random lasing could be achieved 
even in a regime of low scattering. This is a feature of a system exhibiting a high level 
of gain, in which the threshold can be reached even with a low-quality feedback (i.e., 
a cavity with a poor quality factor). 

4-5. Emitted intensity above threshold 

In this section we focus on the characterization of the Mollow random laser above 
threshold. The threshold of a random laser with Mollow gain can be predicted using 
the lasing condition ^ together with the polarizability given by equation ifTTj) . This 
polarizability is obtained in the weak probe Hmit [43]. When the lasing threshold is 
surpassed, Letokhov's theory leads to an exponential growth of the laser intensity 
versus time, and hence steady state cannot be reached. In order to avoid this 
unphysical effect, as already pointed in [21], saturation effects must be included in 
the description of the atomic polarizability at both pump and probe frequencies. 

To do so, optical Bloch equations in the strong probe regime can be numerically 
solved in order to obtain the atomic polarizability. In this case the atomic 
polarizability depends on the lasing intensity in the medium through its associated 

Rabi frequency I^IrlP oc /^l^- Hence, a = a{S, A,n,nYih) (all frequencies will be in 
unit of r in the following equations). The steady-state value of the random laser 
Rabi frequency is considered as the value at which losses exactly compensate gain. 
In this situation, the output power of the random laser -P^™*'' equals the generated 
power in the lasing medium, i.e., P^^^^ oc CTgll^RLp, with dg = ctq (|(5p — Im(Q:)) 
the gain cross-section. On the other hand, the pump-induced fluorescence power is 



Hence, in this regime the ratio of the lasing power to the fluorescence power 



PFiuo(xao|f^|V(l + 4A2 + 2|f]|2). 



Towards a random laser with cold atoms 



19 



induced by the pump can be estimated. This ratio is an important step towards the 
performance characterization of the random laser, as it allows to quantify the amount 
of laser signal that can be extracted from the fluorescence background. From the 
previous results we get 



As can be seen in figure [TOl the ratio of lasing power to fiuorescence background 
is on the order of 5% for a pump detuning A = 1, leading to a measurable signal. 
One interesting characteristic of this random laser concerns the laser frequency at its 
maximum emission power. As shown in figure [TOl this maximum emission frequency 
shifts as the pump intensity increases. This effect is due to the shift in the Mollow 
gain as the pumping power changes. 

4-6. Limitations of the model 

Of course these crude estimations neglect different effects. We assume homogeneous 
pump intensity across the whole system [87] . We are considering one pump frequency 
and one lasing frequency, neglecting thus effects resulting from mode competition |32] 
and inelastic scattering of laser light [92j among others. Despite neglecting those 
effects, these estimations seem to be reasonable at least in recovering the orders of 
magnitude of the lasing threshold and power. 

The inhomogeneities in the atomic cloud density, in the pump intensity and 
in the gain distributions could be incorporated in the RTE model, at the cost of 
a full numerical treatment of coupled equations for the pump beam and emitted 
light intensities. Such an approach has been developed previously in the case of a 
"classical" random laser, i.e., with gain and scattering as separate entities [88]. Another 
improvement would be to deal with a spherical geometry, that is closer to the real 
conditions. Finally, including the Raman gain mechanism into the RTE approach, 
using experimental data as input parameters, would be of great interest in view of 
experiments since the threshold conditions seem easier to reach than with Mollow 
gain, in view of our first estimates. 

5. Outlook and conclusion 

We have presented in this paper our recent investigations on the issue of achieving a 
diffusive random laser in a cloud of cold atoms. We have especially studied different 
gain mechanisms in a situation of standard laser, and we have quantified the random 
laser threshold for two of them. These evaluations show that random laser in cold 
atoms is possible. This is undoubtly our main result. Moreover, we point out that 
the peculiarity of our system would lead to a random laser with a low amount of 
scattering, that is, low feedback. This regime is similar to that encountered in certain 
semiconductor lasers with a very poor cavity, and is different from the working regime 
of random lasers realized to date. 

Moreover, from our experience, we can now draw some preliminary conclusions 
about the comparison of the different gain mechanisms. Application of Letokhov's 
criteria leads to similar critical optical thickness for Mollow gain and Raman gain. 
Nevertheless, the optimum for Mollow gain is obtained outside the range of validity 
of our transport model whereas the amount of scattering is larger with Raman gain. 




(15) 
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Moreover, polarization is better taken into account in our prediction for Raman gain. 
Finally, the optimum pumping parameters are such that pump penetration through 
the cloud will be higher with Raman gain, since detuning and power are larger. For 
all these reasons, Raman gain seems more appropriate to achieve random lasing. 

However, achieving random lasing is not enough, as we need also to detect the 
random laser emission. In our system, gain is almost at the same wavelength than the 
pump (and thus its fluorescence), and this makes the detection challenging. In this 
respect, Raman gain is worse than Mollow gain, as the very small detuning (5r between 
the pump and the gain frequency prevent the use of a Fabry-Perot interferometer to 
distinguish between the laser emission and the pump-induced fluorescence. This could 
be done more easily with Mollow gain (see, e.g., [93]). Moreover, from our preliminary 
evaluations of the intensity of the random laser emission with Mollow gain (section 
14. Sp . we know that the random laser emission should not be too small to be detected. 

For Raman gain we can imagine other ways to extract the optical spectrum of the 
emitted light, by looking at the intensity correlations of the fluorescence, measured 
by a heterodyne technique [94j or a homodyne technique [95j or a correlator [96j. 
Nevertheless these techniques do not enable to filter the random laser light, contrary 
to the use of a Fabry-Perot cavity. Filtering the random laser light could be very useful 
to study the random laser properties, for example the intensity temporal correlations. 
Note that this kind of experiments could be interesting also below the random laser 
threshold [97t[98). 

The third gain mechanism that we have studied is four-wave mixing, which has 
lead to a very efficient cavity-laser. Nevertheless, the configuration we have used so 
far, that is degenerate FWM, is not appropriate at all for the random laser problem, 
as the pump fluorescence and the random laser emission would be exactly at the 
same frequency. However, another conflguration of parametric gain is possible, that is 
non-degenerate FWM, for which two different frequencies and wb are used for the 
two pump fields, and amplification occurs at frequency uj = {cop + (^f) /2. Then, it 
should be possible to choose a large enough detuning to facilitate the detection, and 
to adjust the gain frequency close to the atomic frequency ujo to enhance scattering. 
Preliminary numerical simulations indicate that this configuration is promising [99], 
and its experimental test constitutes our next work. 

Finally, we may also use some additional degrees of freedom that atoms offer. For 
example, cold atomic clouds have been shown to exhibit large Faraday rotations in 
presence of a magnetic field |100) . This Faraday effect may change the threshold of the 
random laser [lOlj . Moreover, magnetic fields can change also the atomic radiation 
pattern [102) and the coherence length of the light transport inside the sample |103j . 
These effects could thus be used to engineer the multiple scattering properties of atoms 
and therefore to achieve a better control on the random laser based on cold atoms. 
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Appendix A. Laser dynamics equations 

We give here the expressions of the functions ^' and 6 that describe the coupUng 
terms in the laser dynamics coupled equations JSHZl) [731174] : 

M 

$(crj, a„) = 2zAi ^ (exp[-i0j](Tja*j - exp[i0j]cr*a„) (A.l) 

m— 1 

M 

+ 2iX2 ^ {exp[-~i(j>j]cr*a'^ - exp[i0j]crja„) + 273 (a^ + a*) , 

m— 1 

M 

'^{aj,a3j,am) = X! (expii^jlAicrsjam + exp[-i(/)j]A2cr3ja;^) (A.2) 

m— 1 

A/ 

- 2iAo ^ (exp[-i0j]o-ja;^ + exp[i(?!)j]crja„) + 70 + 730-3^ - 740-* , 

m— 1 

AT 

e(CT3j, cr-,) = -z^exp[-i0j] (AiCTj - A2(T* + AoCT3j) , (A.3) 
i=i 

where 70 = (r/2)sin26', 73 = (r/4) sin 26* cos 26*, 74 = (r/4)sin2 29, Aq = (.g/4) sin26i, 
Ai = (g/4) (1 + cos 26*) and A2 = (5/4) (1 - cos 261). 

The procedure to numerically solve the equations ([SHZl) follows the method 
presented in [73j. We introduce a family of rt-photon polarizations Sn = 
exp[— in(^j](Tj and S^n = X]j 6xp[~*'^<^j]'''3i- Equations l|5|6p cast then into an 
infinite set of hierarchical equations for Sn and S'3„. In this paper, we truncate this 
set of equations at \n\ = 5, which, in the case of M = 1, leaves us with a total of 
35 equations and corresponds to the generaHzed effective theory GET35, if we use 
terminology of |73j. The hierarchical equations for macroscopic polarizations Sn and 
Ssn, supplemented with equations for amplitudes of modes of the electromagnetic 
field, are then solved numerically. 
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